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ABSTRACT 

We cross-correlate WMAP and ROSAT diffuse X-ray background maps and look for 
common features in both data sets. We use the power spectrum of the product maps 
and the cross-power spectrum to highlight a possible correlation. The power spectrum 
of the product maps does not detect any correlation and the cross-power spectrum 
does not show any significant deviation from 0. Wc explore different explanations for 
this lack of correlation. A universe with a low value of as could naturally explain the 
lack of correlation. We also discuss the systematic effects which can affect this result, in 
particular the subtraction of some cluster signal from the ROSAT diffuse maps which 
could significantly suppress the correlation signal. These systematic effects reduce 
considerably the significance of our constraints on the cosmological model. When we 
include the systematic effects we find a weaker constrain on as allowing models with 
values as large as as — 1 (for fl m — 0.3) to be consistent with the lack of correlation. 
To illustrate the capabilities of the method with future high-quality data, we show how 
from the correlation signal it should be possible to predict the level of contamination 
of the SZ effect on the power spectrum of the CMB. Within the systematic errors 
we find evidence that this contribution is negligible for WMAP and is expected to 
be small in experiments like ACBAR or CBI, but can be important for future high 
resolution experiments. 

Key words: cosmological parameters, galaxies:clusters:general 
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1 INTRODUCTION 

The recent release of the WMAP data (Bennett et al. 2003a) 
has opened a new window for studies of large-scale structure 
based on the well known Sunyaev-ZePdovich effect (SZ ef- 
fect) (Sunyaev & Zel'dovich, 1972). The SZ effect shifts the 
spectrum of the CMB photons to higher frequencies. This 
shift is redshift- independent and proportional to the product 
of the electron column density with the average temperature 
along the line of sight. The electron temperature and opti- 
cal depth due to Thomson scattering are particularly high 
inside galaxy clusters. Thus, the SZ effect is a good tracer 
of clusters, even for those at high redshift. Around galaxy 
clusters, a diffuse, possibly filamentary, distribution of hot 
gas is believed to be present. These filaments have not been 
definitively detected due to their low contrast compared with 
the background (either CMB or X-ray backgrounds). The 
same electrons which cause the SZ effect will also emit X- 
rays (XR) by bremsstrahlung emission. Therefore, one ex- 
pects the SZ effect and the X-ray emission of galaxy clus- 
ters and filaments to be spatially correlated. Since the X- 



ray background and the CMB are not correlated (except at 
very large scales where there could be a correlation due to 
the integrated Sachs- Wolfe effect (Crittenden & Turok 1996, 
Boughn et al. 1998), the cross-correlation of an X-ray map 
with the CMB should enhance the signal of clusters and fil- 
aments with respect to the background. This fact motivates 
the present study. 



We will be interested in studying the cross-correlation 
SZ ® XR (where ® stands for cross-correlation). The main 
advantage of a cross- correlation analysis is that it is possible 
to highlight signals which otherwise would not be detected 
in any of the bands independently (like weak clusters and 
filaments). Therefore, a cross-correlation analysis should, in 
principle, be more sensitive to the low mass range than stan- 
dard studies which first need to detect the objects. This 
may have important implications for studies of the cluster 
physics like pre-heating which are more sensitive to the low 
mass range. Cross-correlation studies can be then considered 
as complementary of the standard analysis based on galaxy 
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clusters. 

To study the cross-correlation signal we need to de- 
fine a statistical object to quantify this correlation. We will 
use the power spectrum of the SZ ® XR map as such 
an object. We will also use the so-called cross-power spec- 
trum (cross-correlation of the Fourier modes). The advan- 
tages/disadvantages of using the power spectrum of the 
product and the cross-power spectrum will be highlighted 
in section 3. 

There are several advantages to using the power spec- 
trum and cross-power spectrum over other statistical ob- 
jects. First, they contain useful information at different 
scales. For instance the mode accounts for the correlation 
coefficient of the two maps. Higher modes will contain infor- 
mation about the fluctuations at smaller scales. The mod- 
elling of the power spectrum is also easier and it can easily 
account for the uncertainties in the assumptions made in the 
model, as we will see below. The power spectrum will also 
tell us something about the contribution of clusters and fil- 
aments to the CMB power spectrum. Previous papers have 
claimed an excess in the CMB power spectrum (Pearson et 
al. 2002; Bond et al. 2002). ft is not yet clear whether this 
excess could be caused by the SZ effect signal or just be inad- 
equately subtracted residuals (compact sources or residual 
noise). An independent estimation of the SZ effect power 
spectrum would help to clarify this point. 

The reader is encouraged to refer to the recent litera- 
ture for a more detailed description of the modelling of the 
power spectrum. In particular he/she may find interesting 
the general discussion given in Cooray & Sheth (2002) , a SZ- 
oriented discussion in Komatsu & Seljak (2002) and Zhang & 
Wu (2003), or an X-ray oriented vision in Diego et al. (2003). 
For the WMAP results, the reader should refer to Bennett 
et al. (2003a) and for ROSAT data he/she can find all the 
relevant information in Snowden et al. (1997). There are also 
several interesting discussions of cross-correlations between 
CMB and X-ray data sets (Kneissl et al. 1997, Boughn et al 
1998), and the expected cross- correlation between WMAP 
and SDSS (Peiris & Spergel 2000). 

She/he may also find interesting the recent results on cross- 
correlating WMAP with radio and X-ray sources looking 
for the 1SW effect (Boughn & Crittenden 2003, Nolta el al. 
2003). Finally, the reader will find interesting the discussion 
in Bennett et al. (2003b) where they find a 2.5a" correlation 
signal between WMAP and the 242 Abell-type clusters in 
the XBACS catalogue (Ebeling et al. 1996). In this work the 
Hubble constant is set equal to 100 h km s _1 Mpc, with h 
generally taken to be 0.7. 



2 WMAP VS ROSAT: CMB VS X-RAYS 

Before starting any description of the model, it is useful to 
give a brief description of the two data sets which are going 
to be used here (the reader should consult the original pa- 
pers for a more detailed description). WMAP data consists 
of 5 all-sky maps at five different frequencies (23 Ghz < v < 
94 Ghz). At low frequencies, these maps show strong galactic 
emission (synchrotron and free- free) . The highest frequency 
maps (41-94 Ghz) are the cleanest in terms of galactic con- 



taminants and will be the most interesting for our purpose. 
The WMAP data is presented in a special format which con- 
serves the size of the pixels and their shape (within small 
deviations) over the sky. This pixelisation (HEALPIX*) is 
very appropriate for power spectrum computation. Within 
this pixelisation, the data is presented with a pixel size of 
w 6.9 arcmin (Nside=512 in HEALPIX ). This pixel size 
over-samples the beam and also is smaller than the pixel 
size of ROSAT. We will repixelise the maps to the next level 
(Nside=256, pixel « 13.75 arcmin). This minimum scale 
(13.75 arcmin) will define a maximum multipole (I = 767) 
beyond which the data does not contain additional informa- 
tion. The units of the WMAP data are temperature fluctu- 
ations with respect to the background (AT). 
We will focus on one basic linear combination of the WMAP 
data, the differenced Q — W bands of the 1° smoothed ver- 
sion of the original data. This differencing completely re- 
moves the main contaminant in this work, the CMB, leav- 
ing a residual dominated by galactic and extragalactic fore- 
grounds as well as filtered instrumental noise. 
On the other hand, the ROSAT All-Sky Survey data (RASS, 
see Snowden et al. 1997) is presented in a set of bands 
(« 0.1 — 2 keV). Low energy bands are highly contaminated 
by local emission (local bubble and Milky Way galaxy) while 
high energy bands show an important contribution from ex- 
tragalactic AGN's. The optimal band for our purposes will 
be the band R6 (f* 0.9-1.3 keV). This band is the best in 
terms of instrumental response, background contamination 
and cluster vs AGN emission. The pixel size is 12 arcmin 
and the units are cts/s/arcmin 2 . The ROSAT maps have 
been cleaned from the most prominent point sources (AGN's 
above 0.02 cts/s in the R6+R7 band). However, we should 
note that for the above threshold (0.02 cts/s), the survey 
source catalogue was not complete and there are still some 
very hard bright sources which were not removed from the 
data. More importantly for our problem is the fact that to- 
gether with the point sources, a fraction of the clusters above 
the 0.02 cts/s threshold were also removed from the maps 
of the diffuse X-ray background. This may introduce a sys- 
tematic error in our interpretation of the signal which we 
should account for. We will come to this point later. 
Due to the different pixel size, we have repixelised the 
ROSAT R6 band using HEALPIX and the same resolution 
level (Nside=256). 

Although the R6 band is the cleanest in terms of galactic and 
AGN contamination, it still contains very strong emission 
coming from the galactic disk. In order to maximise the ex- 
tragalactic signal, we restrict our analysis to regions outside 
the galactic plane. In particular, we will consider only a clean 
portion of the sky above b = 40° and 70° < I < 250° which 
will also exclude the contribution from the north- galactic 
spur. This optimal area of the sky covers w 9% of the sky. 

As mentioned in the introduction, a CMB map will con- 
tain distortions due to the SZ effect and an X-ray map will 
show some structure due to the same hot and dense plasma. 
However, there are many differences between the two emis- 
sion sources which should be well understood before mod- 
elling the power spectrum of the cross-correlation. The dis- 



* available at http://www.eso.org/science/hcalpix. 
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tortions in the CMB map are proportional to the integral of 
the electron density times its temperature along the line of 
sight. When we take the integrated signal across the area of 
the plasma cloud, we find that (assuming T = const), 



Ssz = So 



TM 



(1) 



That is, the total emission depends only on the total pres- 
sure of the plasma cloud, but not on its geometry. The con- 
stant So includes all the proportionality constants (baryon 
fraction, frequency dependence and units, AT/T or mjy). 
On the contrary, the X-ray emission by the same plasma is 
proportional to an integral involving the square of the elec- 
tron density times the square-root of its temperature. If we 
now calculate the total emission from the cloud of plasma 
we find the surprising result that the total emission depends 
very much on the geometry of the cloud. This comes from 
the fact that the bremsstrahlung X-ray emission involves 
two particles and therefore the denser parts of the cloud will 
have a much larger emission rate than the less dense parts. 
Meanwhile, the SZ effect can be very well modelled if we 
only know the amount of gas and its temperature, whereas 
the X-ray emission involves one more unknown degree of 
freedom, the density profile of the electron cloud which is 
poorly known. Actual observations of the X-ray emission in 
galaxy clusters find that the observed total emission cannot 
be simply reconciled with the predictions from analytical 
models. We need to include additional phenomena in the 
models (pre-heating, cooling flows, dumpiness) to explain 
this discrepancy. This suggests that pure modelling of the 
X-ray emission can produce predictions which are far away 
from the observations. In this paper, we will try to overcome 
this problem by modelling the X-ray emission using phe- 
nomenological forms which match the observations. Thus, 
we will model the total X-ray emission as; 



Sxr = 



4ttDi (z 



LoT a (l + z) 
4jr£>,(*) 2 



(2) 



where L x is the X-ray luminosity and the parameters L a ,a 
and tp will be chosen to match the observed L x — T relation. 
In modelling the temperature in both equations (1 and 2) 
we will use the relation, 



T = T M^(l + z)" 



(3) 



The specific values of T ,j3 and <j> will be discussed later. 
The X-ray flux must be converted into the flux units of the 
R6 band. We do this following Diego et al. (2003). 



3 POWER SPECTRUM OF THE PRODUCT 
AND CROSS-POWER SPECTRUM 

The previous discussion relates the mm and the X-ray emis- 
sion from the same plasma. However, our two data sets will 
include other components which could (and eventually will) 
show a spatial correlation between the two maps. A good 
way to highlight this correlation is by using the power spec- 
trum of the product map and compare it with the power 
spectrum of the product of two statistically similar maps 
with no spatial correlation between them. This can be done 
by just rotating one of the maps (so the spatial correlation 
disappears). This approach is different to the standard one 



where one looks for correlations in the Fourier modes (cross- 
power spectrum) . This second approach renders good results 
when the signal responsible for the correlation is extended. 
When one looks for correlations due to compact signals, the 
former approach renders better results. The reason is that 
a cross-correlation of the Fourier modes is equivalent to a 
convolution of the two maps. In this convolution, the spatial 
information of the compact sources is partially lost since it is 
diluted over the Fourier plane. In the absence of noise, both 
approaches should give the same results. However, when the 
noise is present, the correlation between the Fourier modes 
is only evident at large scales since at small scales, their 
correlation produces a signal which is much weaker than 
the oscillations (around 0) of the non-correlated noise. On 
the contrary, by multiplying the two maps in real space we 
make full use of the spatial correlation between the sources 
before going to the Fourier space. We have tested the per- 
formance of the power spectrum of the product against the 
cross-power spectrum of the Fourier modes with simulations 
which try to reproduce the characteristics of our data sets. 
Our results confirm that the power spectrum of the prod- 
uct maps is more sensitive than the standard cross-power 
spectrum. However, the power spectrum of the product has 
one drawback. It is very sensitive to single fluctuations in 
both maps. If we have a 5u fluctuation in each map, then 
after multiplying the final fluctuation will be much larger 
and may dominate the power spectrum. On the contrary, 
the cross-power spectrum is much more stable. In this work 
we will look at both quantities. 

Before modelling the power spectrum of the product 
maps and their cross-power spectrum, it is interesting to dis- 
cuss what else we expect to contribute. We can split our data 
in two components, signal and residual. The signal in our 
case will be the emission (mm or X-ray) of galaxy clusters 
and filaments. The residual will include all the rest. That is, 
the CMB, all the foregrounds, unresolved radio sources and 
the instrumental noise for the case of the WMAP data and 
non-removed AGN's (see above), galactic emission, residu- 
als left after corrections for solar flares, and/or cosmic rays 
plus a small contribution coming from intrinsic instrumental 
read-out noise in the ROSAT case. 

When we cross-correlate the WMAP and ROSAT maps, 
there will be a contribution to the power spectrum com- 
ing from these residuals. Even if the WMAP and ROSAT 
residuals are not correlated, the power spectrum of the prod- 
uct map will show features which are common to some (or 
both) of the residuals. These features will produce a power 
spectrum which will be different from 0. Only the monopole 
(equal to the correlation coefficient) will be if there is no 
correlation between the maps. Then, the fact that the maps 
are not correlated does not mean that the power spectrum 
of the maps must be 0. In order to say whether or not the 
power spectrum of the product map contains a correlated 
signal the easiest way to do it is just to rotate one of the 
maps before we multiply them in real space. If there is a sig- 
nificant correlation between the maps the power spectrum 
of the product of one map times the rotated map will be 
smaller than the power spectrum when the rotation angle is 
0. The power spectrum will reach a minimum for a rotation 
angle larger than the correlation length of the two maps. 
We will call this minimum the background power spectrum. 
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If the maps are not correlated, the power spectrum after 
rotating one of the maps will be roughly the same as the 
power spectrum if we do not rotate them and similar to the 
background power spectrum. 

On the other hand, if the maps are not correlated, the 
cross-power spectrum as a function of the multipole, will os- 
cillate around (with mean value « 0). The identification 
of a correlated signal in the cross-power spectrum is in prin- 
ciple easier since no rotation of one of the maps is required. 
The cross-power spectrum will oscillate around if there is 
no correlation between the maps and will deviate from if 
there is a correlation. 



3.1 Power spectrum of the product maps 

The power spectrum of the product maps is defined as: 



Ct = Cf(2h) + Cf(lh) 



(6) 



'l — < aimCllm > 



(4) 



where the a; m are the spherical harmonics coefficients (and 
their conjugates a ; * m ) of the product map (WMAPxROSAT) 
and the average is over the 2£ + 1 coefficients with the same 
£. This power spectrum can be expressed as the sum of two 
power spectra: 



c 



(5) 



with C® 5 the power spectrum of the product of the resid- 
uals and C? the power spectrum due to the cluster (and 
filament) correlations between the mm and X-ray band. The 
previous equation follows from the assumption that the clus- 
ter and filament signal is not correlated with the residual (all 
the signal which is not due to clusters) . This discussion can 
be illustrated with a simple example. In figure 1 we consider 
a case where the CMB data contains just CMB (simulated) 
and the SZ effect. We cross-correlate this simulated map 
with the real ROSAT R6 data and with a randomised ver- 
sion of ROSAT. The SZ effect emission was simulated based 
on a catalogue of more than 2700 Abell & Zwicky galaxy 
clusters (over « 80 % of the sky). The masses were com- 
puted from the richness and the distances by calibrating the 
magnitude of the 10th brightest member with the known 
distances of 700 clusters. Both masses and distances can 
be very inaccurate when calculated following this process. 
SZ effect total fluxes were computed using equation 1. The 
noise of ROSAT was simulated by randomising the positions 
of the pixels of ROSAT. This technique has the advantage 
that the noise map has exactly the same pdf as the original 
data but no real structure. From figure 1 we can see how 
in fact the cross-correlation of the CMB map plus SZAbM 
with the ROSAT has significant power at small scales. When 
we cross-correlate with the randomised ROSAT this power 
disappears (dotted line). The same thing happens when we 
cross-correlate ROSAT with a 180° rotation of the Abell cat- 
alogue plus CMB (dashed line) . It is important to note that 
even in the case where we cross-correlate CMB + SZAbeii 
with the ROSAT noise, the resulting map still has structure 
at large scales. 

For modelling the term Cf c , we only need to know some- 
thing about the cluster distribution and their signal in each 
band. Basically, this term will be the contribution of two 
terms, 



The first term accounts for the two-halo contribution and it 
includes the contributions to the power spectrum due to the 
cluster-cluster spatial correlation. This term will be signifi- 
cant only at very large scales. However, as we will see later, 
the power spectrum at large scales will be dominated by the 
power spectrum of the cross-correlated residuals, Cf^. Also, 
the large scales will be affected by the window function of 
our optimal area (Sliwa et al. 2001). Therefore, the large 
scales (6 > 20° or £ < 10)) will not be used used here. Since 
the modelling of the two-halo component is a rather compli- 
cated process involving several assumptions about the bias 
and its evolution and that it only contributes significantly 
to the large scales we will not consider the two-halo contri- 
bution in this work. The main contribution at small scales 
will come from the single-halo contribution (C® c (l/i)). This 
is just given by, 



cfAih) 



dM 



dN(M, z) 
dM 



Pi(M,z) 



(7) 



where dV(z)/dz is the volume element, dN(M, z)/dM is the 
cluster mass function and pi(M,z) is the power spectrum 
(multipole decomposition) of the SZ XR cross-correlated 
2D profile of a cluster with mass M at redshift z. In this 
work we will assume Press-Schechter for the mass function 
(Press & Schechter 1974) although other approaches (more 
realistic) could be easily incorporated into the previous for- 
mula. 

The term pi(M, z) can be modelled as 



pi(M,z) =p„(M,z)*f(l,M,z) 



(8) 



where p is just the total signal of the SZ (g) XR cross- 
correlated 2D profile and f(l,M,z) contains the multipole 
dependence which depends only on the geometry of the 2D 
profile. This term can be fitted numerically by the following 
expression (with typical error less than 10 % up to £ = 1000), 



f(l,M,z) = I (exp(-S ; , flc )+cxp(- v ^X)) 



with, 



j2 ^2/(0.97+0. 68e-4/fl c ) 



(9) 



(10) 



where the core radius, R c , is given in rads. Equations (9) 
and (10) are valid for a /3-model with (3 = 2/3 truncated at 
the virial radius. Different /3-models will produce a different 
multipole profile. The effects of the profile will be discussed 
later. The central density is irrelevant for us since we nor- 
malise the total signal using equations 1 and 2. The only 
relevant parameters will be the core radius and the ratio 
p = virial / core radius which we set to p — 10. 
In terms of observable quantities, p can be expressed as, 



p (M,z) = Air\Mean[ 



(11) 



where Mean is the mean signal of the cluster on the sky. 
That is, the sky-averaged product of the mm signal times 
the X-ray signal. 



Ssz(0) = Ssz 



Tot(A) 



Sxr(9) = Sxr 



B(0) 
Tot(B) 



(12) 



(13) 
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Figure 1. The thin solid line is the power spectrum of a cross 
correlated CMB simulation plus a SZ effect simulation (based on 
a catalogue of Abell clusters) with the ROSAT R6 band data. 
In the dashed line we rotate the simulated SZ effect 180 deg in 
the direction E-W. The dotted line is the power of the cross- 
correlation between CMB + SZaoM an d a random realisation of 
the ROSAT data. 



where the terms Ssz and Sxr are given by equations 1 and 
2 respectively. The factors A(9)/Tot(A) and B(6)/Tot{B) 
account for the profile dependence of the signal. It is impor- 
tant to include them because, as compared with the power 
spectrum in the X-rays or the SZ effect (see Diego et al. 
2003), Mean will depend on the assumed profile. From the 
two previous equations, it is easy to show that, 



Mean = 



SszSxr Tot(AB) 
4tt Tot(A)Tot(B) 



(14) 



where Tot(AB) is the integrated 2D profile of the cross- 
correlated SZ (g> XR image while Tot(A) and Tot(B) are 
the integrated profiles of the SZ effect and X-ray 2D pro- 
files respectively. Then, the only additional information we 
need to compute the cluster SZ eg) XR power spectrum is 
to define the scaling relations (equations 1 and 2) and give 
an expression for the core radius as a function of mass and 
redshift. For the scaling relations, we will use the best fit- 
ting model found in Diego et al. (2001). The advantage of 
using this model is that the combinations of parameters of 
this model produce a good fit to several cluster data sets 
(mass function, temperature function, X-ray luminosity and 
flux functions). Later we will discuss other alternatives. For 
the core radius, we will assume that this is given by the 
expression; 



R c 



Rv 
P 



r Mllf i {l + z)- 1 h- 1 Mpc 



(15) 



We will assume that the core radius is a constant fraction 
of the virial radius. We will take this fraction (concentra- 
tion parameter) as p = 10 (Rv = pR c ). We summarise our 
reference model in table 1. We will use this model just for 
illustration purposes. 

Once we have defined our model, we can compute the 
power spectrum (equation 7). 



Table 1. Reference model. All numbers are dimensionless except 
L which is given in units of 10 42 h~ 2 erg /s, T a which is in keV 
and r in h~ 1 Mpc. <j> and tp have been fixed to 1 since we are not 
sensitive to them. This model is in perfect agreement with several 
cluster data sets (Diego et al. 2001, Diego et al. 2003) 





r 


L 


a T 


13 r a 


P 


0.3 0.8 


0.2 


1.12 


3.2 9.48 


0.75 0.13 


10 



3.2 Cross-power spectrum 

The SZ (g> XR cross-power spectrum of (Ce(X)) is defined 
as; 



C t {X) 



, SZ XR* 

<- a im a im 



> 



(16) 



where af^ are the coefficients of the spherical harmonics de- 
composition of the SZ effect map and afj? are the complex 
conjugate of the coefficients of the cluster XR map. The 
modelling of the cross-power spectrum is difficult since it 
involves the direct modelling of the a tm instead of their dis- 
persion, d, however, under certain special conditions this 
complicated modelling can be simplified enormously. 
If we impose that the cluster XR map is proportional to 
the SZ effect map, then their corresponding a ern 's will obey 
the same proportionality and the problem of modelling the 
cross-power spectrum can be solved easily. We have to point 
out that the above situation does not occur in reality but we 
will show how the previous assumption is a good approach. 
If we look at equations 1 and 2 we realise that in the partic- 
ular case where we take T oc M 54 (Nevalainen et al. 2000) 



and L x oc T 



(e.g Markevitch 1998, Arnaud & Evrard 



1999), then, at low redshift the total SZ effect signal is pro- 
portional to the cluster X-Ray flux. If one chooses to use 
a different scaling relation for T — M then the scaling in 
L x — T should be changed accordingly in order to make the 
SZ signal proportional to the XR flux. For instance, a re- 
lation like the standard T oc M 2 ^ would require a relation 
like L x oc T 2 ' 5 which is still a good description of the obser- 
vations for massive clusters. 

When we calculate the flux in the R6 band and transform 
flux to cts/s (see Diego et al. 2003) we introduce an extra de- 
pendence on the cluster temperature which breaks the pro- 
portionality. However, this extra dependence with T is weak 
for clusters above w 3 keV and could be easily compensated 
with a slightly different exponent in the L x — T relation. 
Also, there is a different dependency with redshift in the 
SZ signal and X-ray flux, but at small redshift (z < 0.1, 
D a (z) m Di(z)) the redshift does not play a significant role. 
Also in Diego et al. (2003), the authors shown that the X- 
ray cluster power spectrum is dominated by the low redshift 
population and the intermediate mass clusters (T € [3, 10] 
keV). Under these circumstances, we can make use of the 
above cosmological coincidence (ai m (XR) oc a em (SZ) when 
T oc M 54 and L x oc T 2 ' 85 ) and we can easily model the 
cross-power spectrum. 



a im 



>= KCf z = ^Cf R Cf 



(17) 



Kafm and 



where K is the proportionality constant, a e „ 
we have used the fact that Cf R — K 2 Cf z . The properties 
of the cross-power spectrum are then a combination of the 
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Figure 2. Dependence of the power spectrum with the cosmolog- 
ical parameters (assuming the XR observation is in R6 band and 
the SZ in the Rayleigh- Jeans (R-J) regime). The thick solid line 
is the reference model of table 1. Dashed line shows the change in 
power when we change erg 0.1 units with respect to the reference 
model (<r s = 0.7 bottom, erg = 0.9 top). Thin solid lines show 
the change when we vary fi m 0.1 units (Q m = 0.2 bottom and 
H m = 0.4 top). Dotted lines show the effect of changing T in 0.05 
units, (r = 0.15 top and F = 0.25 bottom). 
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Figure 3. Dependence of the power spectrum with the cluster 
physics (R6 and R-J). The thick solid line is again the reference 
model. The two thin solid lines show the change when we take 
p = 20 (r = 0.065) and p = 5 (r = 0.26)). If we change from 
a = 3.2 to a = 2.7, the power spectrum changes from the solid 
line to the dotted line. If we change fi = 0.75 to fi = 0.56, the 
power spectrum decreases only 20% (not represented). Changing 
L to 0.7L o moves the solid line to the thin long-dashed line. 
Finally, varying T to 0.7T o changes the solid curve to the thick 
sort-dashed curve. 



properties of the individual spectra of the SZ effect and the 
X-ray. These properties have been discussed in the literature 
and we will not repeat them again (Komatsu & Seljak 2002, 
Zhang & Wu 2003, Diego et al. 2003). However, we will 
explore the properties of the power spectrum of the product 
maps in more detail in the next section. 



4 THE POWER SPECTRUM OF CMB®X-RAY 
AS A PROBE 

From the discussion in the previous sections, we have seen 
that we could expect a cluster signal in the the power spec- 
trum of WMAP <g> ROSAT. This signal can be used to 
constrain the cosmological model and/or the cluster physics 
(T — M, L x — T relations, and cluster geometry). Using 
equation 7, we can predict the power spectrum of clus- 
ters for a wide variety of cosmological models and differ- 
ent assumptions about the physics of the plasma. In figure 
2 we show some examples of the dependence of the clus- 
ter power spectrum with the cosmological parameters. The 
dependence with the cluster physics is shown in figure 3. 
The power spectrum shows an important dependence with 
as and Q. m and a weaker dependence with the shape param- 
eter r. This plot illustrates the enormous possibilities of the 
power spectrum as an independent cosmological discrimina- 
tor. The drawback is that the power is also very sensitive to 
the physics of the plasma (figure 3) so one must be very care- 
ful with the selection of the scaling relations and the density 
profile in order to not introduce a bias in the resulting cos- 
mological parameters. However, we can turn this apparent 
problem into a productive way of studying the intra-cluster 
physics. If the cosmological model is known with some ac- 
curacy, then one can use the power spectrum as a way to 
constrain for instance the extension of the plasma cloud. 



From figure 3, it is interesting to see how when the concen- 
tration parameter changes from 5 to 20, the power increases 
a factor 50 (at I ~ 500). This is a unique dependence which 
cannot be observed when one looks at the power spectrum of 
clusters in the mm or X-ray band (Komatsu & Seljak 2002, 
Diego et al. 2003) . Only when we cross correlate these bands, 
we can make evident the dependence of the normalisation 
of the power on the geometry of the cluster (see equation 
14). Also interesting is to see the dependence of the power 
with the scaling relations. In figure 3 we only illustrate the 
dependence with the scaling exponents a and with the nor- 
malisation constants L and T . The dependence with %j> 
and <f> is weak since the power is dominated by low redshift 
clusters. 

It is possible to trace back the dependence of the power 
spectrum on the scaling relations by just looking at equa- 
tions 14, 1 and 2. In the case of L the dependence is 
just Ci oc L 2 a . In the case of T the dependence is a little 
more complicated since it also enters in the band correction 
(B corr = exp(E min (l + z)/kT) - exp(E max (l + z)/kT)) for 
bremsstrahlung, Ci oc (T *B corr ) 2 . The power shows a weak 
dependence with the (3 exponent. A smaller exponent /3 will 
increase the temperature of the clusters with masses below 
Mis = lO 15 ft _1 M0 and will decrease the temperature of 
clusters above that mass. The total luminosity of the clusters 
with M > M15 will also increase as T a . However, this in- 
crease is compensated by the smaller X-ray band-correction 
which peaks at T « 1 keV and decreases for larger tem- 
peratures. The strong dependence of the power with a is 
easier to follow since in this case the temperature does not 
change (and neither does the band-correction). In this case, 
a smaller a will produce a smaller X-ray luminosity (L x ) 
and consequently a smaller power (Ci oc L x ). 
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5 THE WMAP8ROSAT POWER SPECTRUM 
AND CROSS-POWER SPECTRUM 

In order to maximise the signal-to-noise ratio, we have cre- 
ated a new template of WMAP based on a linear combi- 
nation of some of its bands. Since the CMB is frequency- 
independent, two maps at two different frequencies which 
have been filtered with the same beam will contain exactly 
the same amount of CMB per pixel. In our case, the CMB 
is going to be the mayor contaminant so we should try to 
remove it. This can be done easily if we just subtract one 
band from the other. In our case we will subtract the W 
band from the Q band map (both of them smoothed with 
1°). By doing this we will maximise, the SNR of the SZ 
effect with respect to the noisy background. The resulting 
map will have a linear combination of the filtered noise of the 
two bands, plus foregrounds plus the SZ effect. The last one 
will have a signal proportional to the Compton parameter 
times a factor, Correction(Q — W), equal to, 



Correction{Q — W) = 



Avq 



Avw 



(18) 



where the integrals are over the corresponding bandwidths 
(Ais) and f(u) is the well known frequency dependence of 
the SZ effect. On the other hand, since the maps have been 
smoothed, when we cross-correlate the (Q-W) band map 
with ROSAT, we have to keep in mind that, after smooth- 
ing, the 2D profile of the clusters in the CMB map will have 
a 2D profile different from the 2D /3-model. We will also 
include this fact in our calculations. Finally, as we pointed 
out before, we will consider only a clean portion of the sky 
(6 > 40° 70° < £ < 250°) to minimise the correlations in- 
troduced by the galaxy. In this area of the sky we have also 
removed two bright point sources which were not removed in 
ROSAT, MRK 0421 and RBS 0768. In particular we found 
that one of these sources (MRK 0421) was an important 
source of correlation between WMAP and ROSAT (a simi- 
lar correlation was found in Kneissl et al. 1997) . This X-ray 
source is known to be also a powerful radio source. In order 
to avoid being dominated by a single bright source we also 
removed the central part of the cluster A 1367 which shows a 
very intense central emission which could produce spurious 
signals due to fortuitous alignments with the WMAP. How- 
ever, the results proved to be insensitive to the exclusion of 
the central part of this source. 

The power spectrum of (Q-W) WMAP <g) ROSAT is 
shown in figure 4. As we pointed out before, the fact that 
the power spectrum is different from does not mean that 
there is a correlation. To test whether or not the two maps 
are correlated we have to rotate one of them by several rota- 
tion angles and see if there is a trend in the power spectrum 
toward smaller values when we increase the rotation angle. 
After rotating WMAP several degrees (0.2° < 6 < 30°) 
we did not find any significant deviation in the power spec- 
trum of the product map. The power spectrum of the ro- 
tated maps oscillates around 20 % in both directions (up 
and down) for different rotation angles. This means that 
what we observe in figure 4 is just the background power 
spectrum we should see if there is not a significant corre- 
lation between the WMAP and ROSAT maps. We obtain 
the same negative result when we look at the cross-power 



1000 r 




1000 



Figure 4. Power spectrum of (Q-W) WMAP band times ROSAT 
R6. A similar power is observed after rotation of WMAP by sev- 
eral degrees showing that we do not see any significant correlation 
between the maps. Although not shown, the error bars due to map 
rotation are about 20 %. The thin solid lines are the predicted 
signals for the models in table 2. 



Table 2. These three models are compared wit the observed 
power spectrum in figure 4 and the cross-power spectrum in figure 
5. 



Model 




0~8 


L„ 


a 


T 


P 


r a 


P 


A 


0.3 


0.8 


1.5 


2.85 


8.5 


0.54 


0.13 


10 


B 


0.3 


0.9 


1.5 


2.85 


8.5 


0.54 


0.13 


10 


C 


0.3 


1.0 


1.5 


2.85 


8.5 


0.54 


0.13 


10 



spectrum (figure 5). In this case we do not observe any sig- 
nificant deviation from 0. This is not surprising since, the 
cross-power spectrum is less sensitive to correlations at small 
scales. Furthermore, the range of £'s at which is more sensi- 
tive (low £'s) is affected by the window of our selected area 
of the sky. Also, the error bars due to cluster cosmic vari- 
ance are larger at low Cs. In figure 5 we do not show the 
error bars which should include all the above factors and also 
the errors due to the systematic effects (see below). Their 
computation is not a trivial task and they are of no use in 
our case since we can not claim any clear detection in the 
cross-power spectrum. However, it is important to note that 
the error bars in the cross-power spectrum will be less af- 
fected by the systematics than the power spectrum of the 
product maps. For instance, the normalisation of the cross- 
power spectrum will not depend on the assumptions made 
on the internal distribution of the gas (i.e core radius) in the 
clusters as it happens in the case of the power spectrum of 
the product maps. The cross-power spectrum will be then 
an interesting quantity with future data due to its superior 
stability compared with the power spectrum of the product 
map. 

We show a smoothed version of the cross-power spec- 
trum in figure 6. We have rebinned the cross-power spec- 
trum in bins of Al = 20 in order to get a smooth version but 
still is difficult to see any significant deviation from 0. Error 
bars in this case are just the dispersion of the individual 
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Figure 5. Cross-power spectrum (dotted line) compared with the 
three models in table 2. 



data points in each bin. 

The estimation of the cross-power spectrum (and the 
power spectrum of the product maps) at large scales is af- 
fected by large error bars due to the cosmic variance (of 
galaxy clusters) . These error bars are inversely proportional 
to the area of the sky used in the analysis (9 % in our case). 
A small area may contain a cluster population which is not 
representative of the mean cluster population. Cosmic vari- 
ance can be an important problem for the cross-power spec- 
trum but is not a major problem for the power spectrum of 
the maps since it is more sensitive to the small scales where 
our area of the sky is large enough to make the errors due to 
cosmic variance small. Although we do not detect any signal 
neither with the power spectrum of the product maps nor 
the cross-power spectrum, we can still use this fact to set 
some constraints on the model. In figures 4 and 5 we com- 
pare the measured power and cross-power with three differ- 
ent models where we change the parameter as (the models 
are listed in table 2). This simple comparison tell us that 
the only way to accommodate models with high as is by 
reducing the luminosity of the clusters, and/or their tem- 
perature and/or increasing their sizes (decreases the power 
at small scales). They could also be accommodated if the 
SZ effect is significantly contaminated by point sources so 
the net distortion in the CMB is smaller than if the cluster 
signal is just pure SZ effect. 

In the next section we will see how high as can be also 
marginally accommodated if we account for the possible sys- 
tematic effects. 



6 SYSTEM ATICS 

As already mentioned before, one of the main sources of sys- 
tematic error in our conclusions could be the fact that in the 
ROSAT diffuse X-ray background map some of the clusters 
were removed together with the point sources. This could 
contribute to explain the lack of correlation between ROSAT 
and WMAP data. It is difficult to estimate the overall ef- 
fect of this cluster signal subtraction since only a portion of 
the clusters having fluxes larger than 0.02 cts/s (in R6+R7) 



1 .0 
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200 400 600 800 

Figure 6. Binned (A£ = 20) cross-power spectrum of WMAP(Q- 
W) times ROSAT(R6). The cross-power spectrum has been di- 
vided by the factor y/ 'C e (W M AP)C e (ROS AT) to make the fluc- 
tuations at large £ more evident. 
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Figure 7. Systcmatics I. Effect of cluster subtraction in ROSAT. 
The thin solid line corresponds to model C in table 2. The dot- 
ted lines show the reduction in power when we consider that all 
clusters with fluxes larger than 0.02 cts/s (in R6+R7) and sizes 
smaller than 5 arcmin (top), 10 arcmin (middle) and 15 arcmin 
(bottom) are removed from the data. 



were removed from the data. The brightest clusters can still 
be seen in the data. In figure 7 we show the reduction in 
power due to galaxy cluster subtraction for the model with 
as = 1- The top line is the power spectrum when we as- 
sume that no cluster has been removed from the ROSAT 
data (model C). The lines below the top one correspond to 
models where (from top to bottom) we assume that all clus- 
ters with fluxes larger than 0.02 cts/s (in R6+R7) and core 
radius smaller than 5, 10 and 15 arcmin respectively have 
been removed (15 arcmin corresponds roughly to the max- 
imum extent of the PSF of ROSAT). The plot shows how 
we can have overestimated the theoretical power spectrum 
by an order of magnitude at small scales (the cross-power 
spectrum shows a weaker dependence and the reduction in 
power is roughly a factor 3 in the worst case) . 
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Figure 8. Systcmatics II. Effect of the incomplete knowledge of 
the scaling relations. Solid line is the pessimistic case considered 
above (bottom dotted line in figure 7). The dotted lines show 
the effect of pushing down the temperatures, luminosities and 
increasing the size of the clusters. In the top dotted line we reduce 
T to 8.0 keV, in the middle dotted line we also increase the size 
of the clusters to r = 0.16 (and T = 8) and finally in the bottom 
dotted line we also reduce the luminosities to L = 1 (and T = 8, 
r c = 0.16). 



The second source of systematic errors in our conclu- 
sions is in the uncertainties in the cluster scaling relations. 
The models considered in table 2 are just examples of scaling 
relations which more or less reproduce the observations but 
the reality is that the exact form of the underlying scaling re- 
lations is unknown. In figure 8 we try to account for this pos- 
sible source of systematic error by pushing the parameters of 
the models to the limits of the observational constraints. The 
thin solid line corresponds to model C in table 2 and after 
assuming that all clusters with fluxes larger than 0.02 cts/s 
(in R6+R7) and with core radius smaller than 15 arcmin 
have been removed from the ROSAT diffuse maps (bottom 
dotted line in figure 7) . In the dotted lines we change the pa- 
rameters of table 2, T a = 8 (top), T a = 8, r a = 0.16 (middle) 
and T = 8,r = 0.16, L = 1.0 (bottom). The exponents 
a and (3 were already close to the minimum of the observa- 
tional constraints. If we increase them, the power spectrum 
will increase as well (see figure 3). From figure 8 we see that 
models with as ~ 1 are marginally consistent with the lack 
of correlation only when we push down the cluster scaling 
relations to the observational limits and when we consider 
the pessimistic situation where all clusters with fluxes larger 
than 0.02 cts/s (in R6+R7) and with core radius smaller 
than 15 arcmin have been removed from the ROSAT diffuse 
maps. 



7 CONCLUSIONS 

Using the power spectrum of the product maps and the 
cross-power spectrum we do not detect any correlation be- 
tween WMAP (Q-W) map and the diffuse ROSAT X-ray 
background maps. Our estimators (power spectrum of the 
product maps and cross-power spectrum) are dominated by 
the signals of the uncorrelated residuals. Also, the lack of a 



high frequency channel in WMAP does not allow to optimise 
the CMB subtraction while keeping (or even enhancing) the 
SZ signal. The linear combination of Q-W bands removes a 
significant fraction of the SZ signal which makes more diffi- 
cult the detection of any correlation. However, the fact that 
we do not observe a correlation can be used to set limits 
on the model. A low value of as could naturally explain the 
lack of correlation. However, this can be also partially ex- 
plained by the fact that during the point source subtraction 
process in the ROSAT maps (Snowden et al. 1997), a sig- 
nificant fraction of the compact clusters having fluxes larger 
than 0.02 cts/s (in R6+R7) may have been removed together 
with the point sources. This would introduce a bias in the 
theoretical models which would over-predict the correlation 
signal. Therefore, because of the uncertainties in the data 
and the modelling process, the presented constraints are less 
constraining on as than current results from X-ray cluster 
samples. 

Due to the uncertainties in the data and our partial knowl- 
edge of the cluster scaling relations, our constraints should 
be interpreted with caution by a conservative reader. Af- 
ter understanding the different systematics, the reader must 
choose whether or not he/she wants to believe the con- 
straints presented in this work. 

However, to illustrate the capabilities of the method, 
we have considered a pessimistic scenario where the scaling 
relations are pushed down to the observational limits and 
we account for the possible systematic error introduced by 
the cluster removal in ROSAT data we find that high values 
of as (as > 1 for il m = 0.3) seem to be difficult to reconcile 
with the absence of an observed correlation. Due to the sys- 
tematics, our constraints on as are weak but they illustrate 
the capabilities of the method which should render much 
better results with high quality data. 

Since our constraints are basically limited by the un- 
certainties in the data, a similar analysis carried out on the 
4yr WMAP data and on an updated version of the Snowden 
et al. (1997) diffuse X-ray background maps (with a more 
careful point source removal) should render much better con- 
straints or even a detection of the cross-correlation signal. 
Although WMAP's sensitivity to clusters is much poorer 
than ROSAT, the cross-correlation of both maps should en- 
hance weak cluster signals not detected in any of the maps 
individually making the correlation maps an ideal data set 
to look for weak cluster signals. In particular, distant clus- 
ters will produce a weak signal in ROSAT while they can 
produce a significant signal in WMAP. Moreover, an anal- 
ysis based on the cross-correlated maps is more sensitive 
to the cluster model than a similar analysis on the X-ray 
or SZ maps individually. This opens the door to interesting 
studies of the cluster physics using the power spectrum if 
one assumes that the cosmological model is known. In these 
studies, one of the major advantages compared with classi- 
cal studies will be that one can use the raw data since there 
is no need to detect the clusters or to estimate their total 
fluxes extrapolating their observed profile. 

We also have shown how to model the power spec- 
trum of the product map and the cross-power spectrum with 
an intuitive model based on empirical observations (cluster 



© 0000 RAS, MNRAS 000, 1-11 



10 Diego et al. 




10 100 1000 10000 

Figure 9. Current estimates of the CMB power spectrum com- 
pared with predicted SZ effect power spectrum (R-J) for the mod- 
els in table 2. The top solid line is a rebinning (10 bins) of the 
original WMAP CMB power spectrum. The symbols are current 
estimates by CBI and ACBAR (error bars have been omitted ex- 
cept in the last two points) . The last three symbols at I fa 3000 arc 
the estimated power spectrum at high I by CBI (top), ACBAR 
(middle) and the expected CMB power spectrum for a standard 
model (bottom star). Solid lines represent the models in table 2. 
Dotted line has the same parameters as model C but with T = 9 
and r = 0.1. This model should produce a cross-correlation sig- 
nal which is not observed. 



scaling relations) rather than pure modelling of the electron 
density. The cluster power spectrum of CMBcg>X-ray exper- 
iments can be a powerful technique in future cosmological 
studies but can also be useful for studying the physics of the 
intra-cluster plasma. The fact that high values of as seem to 
be difficult to reconcile with the absence of significant cor- 
relation could also be used to rule out the possibility that 
the excess in ACBAR and CBI is due to SZ effect (at least 
for Gaussian models of structure formation). We illustrate 
this point in figure 9 where we compare the power spectrum 
of the SZ effect for the models in table 2 with the recent es- 
timate of the CMB power spectrum by WMAP (solid line) 
and with estimates from ACBAR (Kuo et al. 2003) and CBI 
(Pearson et al. 2002). The model C (ag = 1) does not explain 
the excess in power in the high-£ regime. For comparison we 
also show a model (dotted line) with the same cosmological 
model (as = 1, f2 m = 0.3) but with a higher normalisation 
in the T — M relation (T — 9 keV) and with a smaller 
core radius normalisation (r — 0.1 hT 1 Mpc). This model 
can be compared with the as = 1 models in Bond et al. 
(2002). The dotted line model produces more power than 
model C because of the higher temperature of the clusters 
and their smaller size (which increases the power at small 
scales). However, this model should produce a detectable 
signal in the power spectrum of the product map (similar to 
the top dotted line in figure 8) even in the pessimistic case 
where we consider that all clusters with fluxes larger than 
0.02 cts/s (in R6+R7) and core sizes smaller than 15 arcmin 
are removed from the ROSAT data (and L a = 1). It seems 
difficult to explain the excess in CBI/ ACBAR as due to SZ 
effect with models which does not show a correlation be- 
tween WMAP (Q-W) and ROSAT. At this point it is worth 



mentioning that non-Gaussian models of structure forma- 
tion (e.g Mathis et al. 2003) could make the trick since they 
can produce a low correlation signal in WMAPxROSAT (ba- 
sically given by the low redshift population) while producing 
a high SZ power (also sensitive to the high-redshift cluster 
population) . 
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